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Abstract 
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atomic nuclei. We discuss the status of the nuclear Hamiltonian derived in this framework and 
some recent applications in few-nucleon systems. Nuclear lattice simulations as a new tool to 
address the many-body problem are introduced and some first results based on that method are 
presented. 
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Nuclear physics is on one side an old and well-established science but on the other 
hand a fascinating and new field. This is related to new experimental facilities 
and techniques but even more so to recent developments in theory. Using modern 
high-performance computers, first attempts are made to calculate atomic nuclei 
directly from Quantum Chromodynamics (QCD), the SU(3) co i or gauge theory of 
quarks and gluons [lj. Complementary to this, starting from the ground-breaking 
work of Weinberg [2j, an effective field theory (EFT) approach to the forces be- 
tween two, three and four nucleons has been developed and applied to a variety 
of nuclear bound states and reactions. This EFT is based on the observations 
that (i) nuclei are composed of non-relativistic nucleons (neutrons and protons) 
and virtual mesons and (ii) that the nuclear interactions feature two very dis- 
tinct contributions, long-range one- and two-pion exchanges and shorter-range 
interactions, that can be represented by a tower of multi-nucleon operators. As 
the pion is the pseudo-Goldstone boson of the approximate chiral symmetry of 
QCD (for an introduction, see e.g. Ref. [3]), its interactions with the nucleons 
are of derivative nature and strongly constrained by the available data on pion- 
nucleon scattering and other fundamental processes. However, in harmony with 
the principles underlying EFT (for an introduction, see e.g. Ref. [4]), one has 
also to consider operators of nucleon fields only. In a meson-exchange model of 
the nuclear forces, these can be pictured by the exchanges of heavier mesons like 
<7, p, lj, and so on - but such a modeling is no longer necessary and also does not 
automatically generate all structures consistent with the underlying symmetries. 
Also, in the EFT approach, the forces between three and four nucleons are gen- 
erated consistently with the dominant two-nucleon forces - which could never be 
achieved in earlier modeling of these forces. 

Due to the non-relativistic nature of nuclei, the underlying equation for the 
nuclear ^4-body system (where A is the atomic number) to be solved is the 
Schrodinger equation, where the various contributions to the nuclear potential 
are organized according to the power counting discussed below. In a second step, 
bound and scattering states are calculated as solutions of this equation. This 
does not only allow to pin down the various low-energy constants related to the 
multi-nucleon interactions, but also to check the convergence of the approach by 
including higher orders in the underlying potentials. Once data for the two- and 
three-nucleon systems are described with sufficient precision, one is then in the 
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position to perform ab initio calculations of nuclei, eventually combining well 
developed many-body techniques with the forces from chiral nuclear EFT. An- 
other venue to approach light and medium-heavy nuclei is based on simulation 
techniques, as they are so successfully utilized in lattice QCD [5] to calculate the 
properties of protons, neutrons and many other hadrons. All this is accompanied 
and extended by the construction of the corresponding electro-weak charge and 
current operators that allow for many further fine tests of the structure of nuclei 
and also of the calculation of fundamental nuclear reactions that are of relevance 
to the generation of the elements in the Big Bang and in stars. 

Disclaimer: Clearly, this is not a detailed all purpose review of this field but it 
rather intends to give an introduction to the underlying ideas and some recent 
applications. Two recent detailed review articles are Refs. [6j|7j that contain also 
many references to earlier work. We have therefore not attempted to be complete 
or exhaustive in the references, but the interested reader will find sufficient quotes 
to the literature to be able to acquire a much deeper understanding. 

Our article is organized as follows: Section 2 contains the basic ideas of the 
chiral Lagrangian and power counting of the nuclear forces. We end this section 
with a brief review of the current status of the nuclear Hamiltonian derived in this 
framework. Section 3 contains some applications of these forces to nuclei, based 
on calculations using exact few-nucleon methods. In section 4, the new method 
of nuclear lattice simulations (nuclear lattice EFT) is presented and some first 
results obtained in that scheme are displayed. 



2 From the effective chiral Lagrangian to nuclear forces 

Our goal is to develop a systematic and model-independent theoretical frame- 
work capable to describe reactions involving several nucleons up to center-of-mass 
three- momenta of (at least) the order of the pion mass M n . Following the usual 
philosophy of effective field theory, we aim at the most general parameterization of 
the amplitude consistent with the fundamental principles such as Lorentz invari- 
ance, cluster separability and analyticity. Given that the energies of the nucleons 
we are interested in are well below the nucleon mass, it is natural and appropriate 
to make use of the non-relativistic expansion (i.e. an expansion in inverse powers 
of the nucleon mass mjv). Accordingly, in the absence of external probes and 
below the pion production threshold, we are left with a potential theory in the 
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Figure 1: Singularity structure of the partial- wave two-nucleon scattering am- 
plitude in the complex energy plane. The solid dot indicates the position of 
the S-wave (virtual) bound state. The shaded area shows the region where the 
effective range function k 2l+1 cot^(fc) is a meromorphic function of k 2 . 

framework of the quantum-mechanical j4-body Schrodinger equation 



The main task then reduces to the determination of the nuclear Hamilton oper- 
ator Hq + V. This can be accomplished using the framework of chiral perturba- 
tion theory (ChPT) [2]. Notice that the approach outlined above automatically 
maintains unitarity of the scattering amplitude and correctly reproduces its an- 
alytic properties at very low energies. Consider, for example, the singularities of 
the neutron-proton / th partial wave amplitude in the complex energy plane with 
E = k 2 /mjV) visualized schematically in Fig. [I] The discontinuity across the 
right-hand cut running from E = to E = +00 and separating the physical and 
unphysical sheets is determined by elastic unitarity which is already built in the 
Lippmann-Schwinger equation. On the other hand, the left-hand cuts are gov- 
erned by the properties of the interactions. Exploiting only the knowledge of the 
finite-range nature of the nuclear force and parameterizing the nucleon-nucleon 
(NN) potential by zero-range terms is sufficient to correctly describe the analytic 
structure of the amplitude within the region near threshold limited by the branch 
point on the first left-hand cut associated with one-pion exchange and located 
at k 2 = M 2 /A. The so-called pionless EFT, see [8j for a recent overview, is 
based on the Lagrangian involving all possible zero-range two- and more-nucleon 
operators with increasing number of derivatives and can, in the two-nucleon sec- 
tor, be matched to the well-known effective-range expansion, i.e. an expansion of 
the function k 2l+1 cot 5i{k) in powers of k 2 . Clearly, this approach is limited to 
very low energies corresponding to nucleon momenta well below the pion mass, 
see Fig. [TJ To extend the range of applicability, one needs to correctly describe 
the left-hand singularities of the amplitude and thus to explicitly include in the 
potential the contributions emerging from the exchange of one or several pions. 





(1) 
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We will describe in the next subsections how this task can be accomplished in 
a systematic way by exploiting the spontaneously broken, approximate chiral 
symmetry of QCD. Another method to construct the NN scattering amplitude is 
dispersion theory, see Ref. [9j for a recent application of this approach. 



2.1 Chiral Lagrangian and power counting 

Within the framework of ChPT, nuclear forces are derived from the most gen- 
eral effective chiral Lagrangian by making an expansion in powers of the small 
parameter q defined a^j] 

where Q ~ \k\ ~ M n is a typical external momentum (the soft scale) and A 
is a hard scale. Appropriate powers of the inverse of this scale determine the 
size of the renormalized low-energy constants (LECs) in the effective Lagrangian. 
Notice that once renormalization of loop contributions is carried out and the 
renormalization scale is set to be fi ~ M v as appropriate in ChPT, all momenta 
flowing through diagrams appear to be, effectively, of the order ~ M n |10|. Con- 
sequently, one can use naive dimensional analysis to estimate the importance of 
(renormalized) contributions of individual diagrams. 

To be specific, consider a connected Feynman graph with iV nucleon linesj^] It 
is easier to count the powers of the hard scale A rather then of the soft scale Q by 
observing that the only way for A to emerge is through the corresponding LECs. 
Thus, the low-momentum dimension v of a given diagram can be expressed in 
terms of the canonical field dimensions Ki + 4 of Vi vertices of type i via 

v = -2 + s ^V i K i , ^ = di + -rii+pi - 4, (3) 

where rii (Pi) an d di refer to the number of the nucleon (pion) field operators 
and derivatives or pion mass insertions, respectively. The constant —2 in the 
expression for v is just a convention. The power counting can also be re-written 
in terms of topological variables such as the number of loops L and nucleon 
lines N rather than Ki which are appropriate for diagrammatic approaches. For 
connected diagrams the above equation then takes the form 

v = -A + 2N + 21 + ^^^, Ai = di + -ni-2. (4) 

Note that we use this parameter and the soft scale Q synonymously in what follows. 
2 We remind the reader that nucleons cannot be destroyed or created within the nonrelativistic 
approach. 
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Chiral symmetry of QCD guarantees that Goldstone Bosons which, in the case of 
two light flavors are identified with the pions, couple only through vertices involv- 
ing derivatives or powers of M n . This implies that the effective Lagrangian con- 
tains only irrelevant (i.e. non-renormalizable) interactions with m > 1 (Aj > 0) 
which allows for a perturbative description of pion-pion and pion-nucleon scat- 
tering as well as nuclear forces. The leading interactions, i.e. the ones with the 
smallest possible Aj, that is Aj = 0, have the form 



£<°> = -d ll 7T-d^7V--M^ 2 + N^ 

2 2 71 



N 



~ ^C s (N j N)(N j N) - -Cx(N^aN) ■ (N j aN) + ... , (5) 

where 7r and N refer to the pion and nucleon field operators, respectively, and a 
(t) denote the spin (isospin) Pauli matrices. Further, g& (F„) is the nucleon axial 
coupling (pion decay) constant and Cs,t are the LECs accompanying the leading 
contact operators. The ellipses refer to terms involving more pion fields. It is 
important to emphasize that chiral symmetry leads to highly nontrivial relations 
between the various coupling constants. For example, the strengths of all Aj = 0- 
vertices without nucleons with 2, 4, 6, . . . pion field operators are given in terms 
of F-x and M n . Similarly, all single-nucleon Aj = 0-vertices with 1,2,3,... pion 
fields are expressed in terms of just two LECs, namely qa and F n . We refer 
the reader to Refs. [11[|12| for further details on the construction of the effective 
chiral Lagrangians (for a modern way to construct the pertinent pion-nucleon 
Lagrangian, see e.g Ref. |13|). 

The expressions for the power counting given above are derived under the 
assumption that there are no infrared divergences. This assumption is violated for 
a certain class of diagrams involving two and more nucleons due to the appearance 
of pinch singularities of the kind 

to + ie Lq — ic 

Here, i/(lo + ie) is the free nucleon propagator in the heavy-baryon approach cor- 
responding to the Lagrangian in Eq. ^ . Clearly, the divergence is not "real" but 
just an artefact of the extreme non-relativistic approximation for the propagator 
which is not applicable in that case. Keeping the first correction beyond the static 
limit, the nucleon propagator takes the form i/(lo — I 2 / (2mjv) + ie) _1 leading to 
a finite result for the integral in Eq. ([6]) which is, however, enhanced by a fac- 
tor ni]\f/\q\ as compared to the estimation based on naive dimensional analysis. 
In physical terms, the origin of this enhancement is related to the two-nucleon 
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Green's function of the Schrodinger equation 0. The nuclear potential V we 
are actually interested in is, of course, well defined in the static limit — > oo 
and thus not affected by the above mentioned infrared enhancement. The precise 
relation between the nuclear potentials and the amplitude corresponding to a 
given Feynman diagram will be discussed in the next section. 

It is now instructive to address the qualitative implications of the power count- 
ing in Eq. Q and the explicit form of the effective chiral Lagrangian. First, one 
observes that the dominant contribution to the nuclear force arises from two- 
nucleon tree-level diagrams with the lowest-order vertices. This implies that the 
nuclear force is dominated by the one-pion exchange potential and the two contact 
interactions without derivatives. Pion loops are suppressed by two powers of the 
soft scale. Also vertices with A, > involving more derivatives are suppressed 
and do not contribute at lowest order. One also observes the suppression of many- 
body forces: according to Eq. Q, iV-nucleon forces start contributing at order 
Q~ 4+2N . This implies the dominance of the two-nucleon force with three- and 
four-nucleon forces appearing as corrections at orders Q 2 and Q 4 , respectively. 



2.2 Derivation of the nuclear forces 

We now clarify the meaning of the nuclear potential V and outline some ap- 
proaches that can be used to derive it. We first rewrite Eq. ([I]) as the Lippmann- 
Schwinger (LS) equation for the half-shell T-matrix 

T af 3 = Vafi+^2 F i ; t T 7^ ( 7 ) 

^ tip- b^ + ie 

where a, (3 and 7 denote the few-nucleon states, E a is the kinetic energy of the 
nucleons in the state a and ^ is to be understood as a sum (integral) over all 
discrete (continuous) quantum numbers of the nucleons. On the other hand, the 
scattering matrix S, which is related to the T-matrix via 

S af5 = 5(a -p)- 2niS(E a - E p )T aP , (8) 

can be directly computed from the effective chiral Lagrangian using the Feynman 
graph technique. This then allows one to define the potential V Q p by matching 
the amplitude to the iterative solution of Eq. ([7]) which in the operator form can 
be written as 

f = V + VG V + VG VG V + ... , (9) 

where Go is the A-nucleon resolvent operator. The outlined approach is, of course, 
not new and has been extensively used in the fifties of the last century in the 
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context of the meson field theory, see e.g. 14,15]. There is, however, one subtlety 
here related to the fact that the T-matrix calculated from the effective Lagrangian 
using Feynman diagrams and entering Eq. Q is only available on the energy shell. 
On the other hand, the potential to be substituted in the LS equation is needed 
off the energy shell. This results in an ambiguity in the definition of the potential 
corresponding to the freedom in carrying out off-the-energy-shell extension. This 
should not come as a surprise given that the Hamiltonian Hq + V is not an 
observable quantity. 

An alternative method to define nuclear forces exploits another old idea of 
decoupling the pion states from the rest of the Fock space by means of a suitably 
chosen unitary transformation [16|. This approach was formulated in the context 
of chiral EFT in Refs. (TTj. The derivation of the unitary operator, nuclear forces 
and currents can be carried out straightforwardly using perturbation theory in 
powers of Q and employing the "algebraic" version of the power counting in 
Eq. ([3]). The above mentioned ambiguity of the nuclear potentials and currents 
can be systematically explored in this approach by performing further unitary 
transformations after decoupling the pion states. Interestingly, this ambiguity 
turns out to be strongly constrained by renormalizability of the Hamiltonian |18| . 

To be specific, consider the derivation of the long-range two-nucleon potentials 
up to next-to-next-to- leading order (N 2 LO). For the sake of simplicity, we use 



here the matching approach and follow closely Ref. 19 , see also Ref. 20 for a 
pioneering calculation within the framework of time-ordered perturbation theory. 
We will not consider here the short-range part of the nuclear force since it can 
be directly read off from the Lagrangian. Notice that it is not necessary to 
explicitly evaluate pion loop diagrams involving contact interactions unless one 
is interested in the quark-mass dependence of the short-range operators. As 
already pointed out before, the leading-order (LO) contribution ~ Q° is due to 
one-pion (l7r) exchange. Evaluating the contribution of diagram (a) in Fig. [2] for 
on-shell nucleons yields 

2 



t/(0) _ T (o) 

V 1tt — 1 1tt 



9 A \ o\ ■ qo 2 ■ q , . 



where q = p' —pis the nucleoli momentum transfer, while p, p' refer to the center- 
of-mass (CMS) initial and final momenta. Notice that while V-f® is uniquely 
defined in this procedure in the static limit with Ep = Ep/ = 0, the relativistic 
corrections are not since one is, in principle, free to add to V\ n terms proportional 
to — Epi . 
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(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



(h) 



Figure 2: Diagrams contributing to the long-range part of the two-nucleon poten- 
tial at LO (a), NLO (b)-(e) and N 2 LO (f), (h). Solid and dashed lines represent 
nucleons and pions, respectively. Solid dots and filled rectangles refer to the 
leading (Aj = 0) and subleading (Aj = 1) vertices from the chiral Lagrangian. 



Because of parity conservation, next-to-leading order (NLO) corrections to the 
potential appear at order Q 2 rather than Q \2\. For all two-pion (2tt) exchange 
diagrams in Fig. [2] except the box graph (e), the potential can be defined via the 



identification K 



(2) 



2tt 



(2) I 

ToJ ^ ^ . In the case of the box diagram, we have to 



subtract the iterated one-pion exchange contribution V^ 1 GqV^ to avoid dou- 
ble counting. Evaluating the corresponding Feynman diagram in the CMS, one 
obtains the contribution proportional to the integral 

dH (2m N ) 2 i 



(2vr)4 [{p-lf 
dH 



N 



N 



+ ie] [l\ — M 2 + ie] [I 2 — M 2 + ie) 



1 



(2tt)3 \^(Ep 



E p-T+ ie ) w 2 



2ujlUJ%((jJi + UJ2) J 



where oji = y I 2 + M 2 and the virtual pion momenta are given by l\ =l,h = l+Q- 
Notice that we used here the relativistic expressions for the nucleon propagators 
in order to avoid the pinch singularity discussed in the previous section. The 
first term in the round brackets is nothing but the iterated one-pion exchange, 
V-f£ GoV\n , while the second one gives rise to bax . Notice further that, as 
explained above, the iterated lir exchange is enhanced compared to the estimation 
~ Q 2 based on naive dimensional analysis since E ~ 0(Q 2 /m^) <C 0(Q). 

The final result for the 27r-exchange potential at order Q 2 can now be ob- 
tained by evaluating the loop integrals and carrying out the spin-isospin algebra. 
Employing dimensional regularization and using the decomposition 

V = Vc+T\-T2 Wc+[Vs + Ti • r 2 Ws] di-a2+[Vr + ti • r 2 Wt) arq°2-q, (12) 
it takes the form 

(2) _ M\q\ 



c 

(2) 



384tt 2 F4 

iyP) - 
q 2 s ~ 



AM 2 (hg\ - Ag\ - 1) + q\2Zg\ - 10 5 1 - 1) + 



3^1 



64tt 2 FJ 



L(W\), 



4M 2 + q 2 

(13) 
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where we do not show the contributions that are polynomial in momenta as these 
can be absorbed into the contact interactions. The loop function L is defined via 



^4M 2 + q 2 ^ y/4M% + q 2 + \q\ 



(14) 



\q\ 2M^ 
Notice that the ultraviolet divergences entering the loop integrals are polynomial 
in the external momenta and, therefore, do not affect the non-polynomial pieces 
if one uses dimensional regularization or equivalent schemes. 

The next-to- next-to-leading order (N 2 LO) corrections emerge at order Q 3 . 
Again, parity conservation forbids any l7r-exchange contributions at this order. 
The 2-7r-exchange terms emerge from triangle diagram (f) in Fig. [2] through the 



identification V n 



(3) 



2- 



T. 



(3)| 



finds 



p p' 



while graph (h) yields a vanishing result. One 



V, 



(3) 



16yrFi 



it; 



(3) 



1 



V2 W S 



(3) 



2M 2 (2ci - c 3 ) - c 3 q 

,2 



{2Ml + q 2 )A{\q\), 



-^_c 4 (4M^ + q")A(\q\), 



where the c, are LECs associated with the ttttNN vertices of order Aj 
the loop function A is given by 



A{\q\) 



1 



arctan 



\q 



(15) 
1 and 

(16) 



2\q\ 2M W 

In addition to the static terms, there are, in principle, also 1/mjv-corrections to 



NLO graphs, see e.g. the last term in the brackets in Eq. ( 11 ). The nucleon mass 
is, however, often treated as a very heavy scale with ^> A, see Ref. [2] for a 
discussion, leading to a suppression of the 1/mjv corrections. 

It is now instructive to address the convergence of the chiral expansion for the 
long-range two-nucleon force. Given that the obtained expressions depend solely 
on the momentum transfer q, the potential is expected to be local in coordinate 
space: 



V(f) = Vc+ti-t 2 W c + 



Vs + T\ ■ T 2 Ws 



(TV (72 + 



Vt + T\ ■ T2 Wt 



(17) 



Sl2 

where r is the distance between the nucleons, S12 = (3ai -f 02 -r — o\ -^r 2 ) /r 2 and 
Vx, Wx are scalar functions of r = \ f\. The Fourier transform of the expressions 



in Eqs. (13) and (15) is, however, ill-defined as the potentials are not bounded 

\q\. At finite distances r > 0, the potential can be 



as q increases, where q 
obtained through a suitable regularization 

d 3 q 



V{f) 



lim 

A— >oo 



I 



(2tt) s 



e-^V(q)F K {\q\) , 



(18) 
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Figure 3: Chiral expansion of the isovector-tensor (upper row) and isoscalar 
central (lower row) long-range potentials Wt(t) and Vc(r), respectively. The 
left (right) panel shows the results for the EFT without (with) explicit A(1232) 
degrees of freedom. The light-shaded band shows the estimation of the intrinsic 
model dependence associated with the short-range components as explained in 
the text (only shown for the theory without deltas). 



where the regulator function F\(x) can e.g. be chosen as F^(x) = exp(— x 2 /A 2 ). 
Alternatively and more elegantly, one can write the functions Wx and Vx in terms 
of a continuous superposition of Yukawa functions which can easily be Fourier 



transformed, see Ref. 19 for more details. For example, for central potentials 
one obtains the unsubtracted dispersive representation 

2 f°° . ftrdi) ... J 



V C (q) 



dfifi-^ 



V c (r) 



dniie-vrpcfa), (19) 



where pc(fJ-) = I m [Xc(0 + — ip)] is the corresponding spectral function. 

In Fig.|3]we show the chiral expansion for the two most important cases, namely 
for the isovector-tensor and isoscalar-central potentials Wt{t) and Vc{r). We 
also include the contributions at next-to-next-to-next-to-leading order (N 3 LO) 



whose explicit form can be found in Ref. 21 but restrict ourselves to the local 
pieces omitting the 1/mjsr corrections. The shaded bands in the figure visualize 
the estimated scheme dependence which is intrinsic to the separation between 
the long- and short-range contributions in the potential. Specifically, we only 



include in the dispersive integrals in Eq. ( 19 ) the components in the spectrum with 



p < A = 1 GeV. The high-p components generate terms which, at low momenta, 
are indistinguishable from contact interactions parameterizing the short-range 
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part of the chiral potential. The bands correspond to the variation of A in the 
range from 800 MeV to oo. Their widths may, therefore, serve as an estimation of 
the size of short-range components not associated with the dynamics of Goldstone 
Bosons. 

The potential in the isovector-tensor channel is clearly dominated by one-pion 
exchange V\ % . Two-pion exchange contributions in this channel become visible 
at distances of the order r ~ 2 f m and smaller. The strong, attractive isoscalar- 
central potential of intermediate range is another well-known feature of the two- 
nucleon force. Phenomenologically, it is attributed to the correlated two-pion 
exchange which is often modeled in terms of the cr-meson exchange [22 J. In chiral 
EFT, on the other hand, all low-energy manifestations of the a and other heavy 
mesons are systematically taken into account through values of the LECs in the 
effective Lagrangian. The resulting strength of Vq turns out to be comparable to 
that of Vx-k even at distances r ~ 2 fm and appears to be an order of magnitude 
bigger than the strength of two-pion exchange in any other channel. The large 
size of the N 2 LO contributions can be traced back to the large numerical factor 
of 3/16 in the case of Vq \ an enhancement by one power of 7r, cf. Eqs. (13) and 



(15), and the large values of the LECs 03,4. For example, the determination from 
pion-nucleon subthreshold coefficients at order Q 2 leads to C3 = —3.9 GeV -1 



and c 4 = 2.9 GeV" 1 23,24 which are larger in magnitude than the expected 
natural size 1 03^4 1 ~ <?a/A ~ 1 GeV -1 . Moreover, even larger values are obtained 
from pion-nucleon scattering at order Q 3 where effects of pion loops are taken 
into account 25 , 26 . The large values of 03,4 can be traced back to the implicit 



treatment of the A(1232) isobar. Given the fairly low excitation energy of the 
Delta toa — mjv, which is numerically ~ 2M n , one may expect that its explicit 



inclusion in the EFT within the so-called small scale expansion 27 based on the 
(phenomenological) extension of the counting in Eq. ^ to 

\ A ' A ' A j ' l2Uj 

will allow one to resum a certain class of important contributions leading to a 
superior convergence as compared to the delta-less theory. The improved con- 



vergence is indeed observed both in the case of pion-nucleon scattering 28 and 



nuclear forces [24,29,30 . In particular, the dominant contribution to Vc and 
Wt emerges in the A-full theory already at NLO with the N 2 LO contributions 
providing fairly small corrections, see the right panel in Fig. [3} For example, the 
single A-excitation in the box diagrams (d) and (e) of Fig. [2] generates at NLO 
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the isoscalar-central potential 



where Ha denotes the irNA axial coupling. In the standard A-less approach based 
on the assignment ttia — m^v ~ A S> M n , this numerically large contribution 
is shifted to N 2 LO where it is reproduced through the A-isobar saturation of 



C3, = — 4:h A /(9(m^ — m^)) 23 , see Eq. (15). Having included the effects 
of the A-isobar explicitly, one finds strongly reduced values of the LECs 03,4 in 
agreement with the naturalness assumption. For example, using = 3g y 4/(2v / 2) 
from SU(4) or large N c , one obtains c 3 = -0.8 GeV" 1 and c 4 = 1.3 GeV" 1 §5]. 
Thus, the major part of the unnaturally large subleading 2-7r-exchange potential 
at N 2 LO is shifted to NLO in the A-full theory. This more natural convergence 
pattern is visualized in the right panel of Fig. [3j 

When substituted into the Schrodinger equation, the long-range potentials in- 
troduced above provide an approximate representation of the nearby left-hand 
singularities in the partial wave amplitude as shown in Fig. [T] These cause a 
rapid energy dependence and imply non-trivial relations between the coefficients 
in the effective-range expansion which can be regarded as low-energy theorems 



(LETs) [31-33 and confronted with the data. A pedagogical introduction to the 
LETs and their relation to the so-called modified effective range expansion can 
be found in Ref. |34| . Assuming perturbativeness of the pion-exchange contribu- 
tions, these relations can be worked out analytically within the scheme proposed 
by Kaplan, Savage and Wise in Ref. (35J. The resulting LETs, however, appear 
to be strongly violated in the 1 So and the 3 Si- 3 -Di channels 31 32 . This in- 
dicates the non-perturbative nature of the one-pion exchange potential at least 



in these channels, see also Ref. 36 . Employing a non-perturbative treatment of 



the pion-exchange potential, the LETs were tested numerically in Ref. 37 . A 



closely related approach is followed by Birse et al. in Refs. [38 -41 by analyzing 
the energy dependence of the residual short-range potential in a given partial 
wave. Perhaps the most impressive evidence of the chiral 27r-exchange comes 
from the Nijmegen partial wave analysis of proton-proton scattering [42] . Here, 
the Schrodinger equation is solved for a specific choice of the long-range poten- 
tial outside of some boundary b. Short-range physics is then taken into account 
by choosing appropriate boundary conditions at r = b. The number of parame- 
ters needed to describe experimental data below the pion production threshold 
with Xdatum ~ 1 can be- regarded as a measure of the amount of physics not 
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included in the long-range potential. With b = 1.4 fm, the authors of Ref. 42 



observe a reduction 31 — > 28 — > 23 in the number of parameters when employing 

(2) (2) (3) 

Vitt — > V\-k + ~^ Vin + ^271- + ^27r as the long-range potential (in addition to 
the corresponding electromagnetic interactions). 



2.3 Nuclear forces: the status and open issues 

We now summarize the current status of the nuclear forces within the heavy- 
baryon, A-less formulation based on the power counting of Eq. Q. In this 
scheme, the nuclear Hamiltonian is presently worked out up to N 3 LO in the 
chiral expansion, 

H = H + V 2N + F 3 N + ^4N + • • • (22) 



with 



2N ^ V 2N ' V 2N ' V 2N 



V2N 

^3N - V 3N + V 3N + . . . , 

V4N = ^ 4 ( N } + • • • . (23) 

where the ellipses refer to terms beyond N 3 LO. For two nucleons, it turns out 
to be necessary and sufficient to go to N 3 LO in order to accurately describe 
the neutron-proton and proton-proton phase shifts up to laboratory energies of 
Eiab ~ 200 MeV (37,43 . This is visualized in Fig. [4] where, as a representative 
example, the experimental data for the neutron-proton differential cross section 
and vector analyzing power at = 50 MeV are shown in comparison with the 
calculations based on the chiral NN potentials of Refs. |37|[43 and various modern 
phenomenological potentials. At the accuracy of N 3 LO, it is mandatory to take 
into account isospin-breaking (IB) contributions. The dominant IB effects emerge 



from the charge-to-neutral pion mass difference in the l7r and 2ir exchange 47 



(charge-independence breaking), proton-to- neutron mass difference in the 2tt- 
exchange |48[|49) (charge-symmetry breaking) and the two derivative-less contact 
interactions (both charge-independence and charge-symmetry breaking in the two 
S- waves). The short-range part of the potential at N 3 LO receives contributions 
from 24 isospin-invariant and 2 IB contact interactions whose strength was ad- 
justed to phase shifts (scattering data) in Ref. [37] (Ref. [43])- Both available 
versions of the N 3 LO potential employ a finite momentum-space cutoff in order 
to regularize the Schrodinger equation. This cutoff is varied in Ref. [37] in the 
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Figure 4: Neutron-proton differential cross section (left panel) and analyzing 
power (right panel) at E\ a b = 50 MeV calculated using chiral EFT, the CD Bonn 
2000 potential of Ref. 44] and the potential developed by Gross and Stadler in 
Ref. |45| . Also shown are results from the Nijmegen partial wave analysis |46| . 
References to data can be found in 1461. 



range A = 450 . . . 600 MeV. More details on the construction of chiral potentials 
at N 3 LO can be found in the comprehensive review articles [7 



50 



While the chiral expansion of the long-range nuclear forces emerges rather 
straightforwardly, the power counting for the short-range operators and the closely 
related issue of non-perturbative renormalization of the Schrodinger (or LS) equa- 



tion are still being debated in the community, see e.g. 51-53 . Here, the main 



conceptual difficulty is associated with the non-perturbative treatment of the Itt- 



exchange potential V]_ w . In spite of a considerable effort, see e.g. Refs. 35,54 



no approximation to is presently known which would capture the relevant 
non-perturbative physics and, at the same time, be analytically resummable 
and renormalizable. One is, therefore, left with the numerical solution of the 
Schrodinger equation for appropriately regularized chiral potentials along the 
lines of Ref. [55] . This paper also explains the meaning of renormalization in 
such an approach and provides a tool to verify its consistency a posteriori by 
means of the so-called Lepage-plots. Notice that iterating the potential, trun- 
cated at a given order of the chiral expansion, in the LS equation necessar- 
ily generates higher-order contributions in the amplitude, which are, generally, 
ultraviolet-divergent and whose renormalization requires counter terms beyond 
the truncated potential. It is, therefore, not legitimate in such an approach to 
arbitrarily increase the cutoff A. This point is exemplified in Ref. 33 using an 



exactly solvable analytical model, see also Ref. [55] for a qualitative discussion. 
More work is needed in order to (better) understand the power counting for the 
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(a) 



(c) 




(f) 



Figure 5: Various topologies contributing to the 3NF up to order Q 5 . Shaded 
blobs represent the corresponding amplitudes. 



scattering amplitude in the presence of the long-range pion-exchange potentials. 
A promising tool to address this question is provided by the modified effective 



range expansion, see the discussion in Ref. 56 and Refs. 39 40] for related work. 



Also, it remains to be seen whether renormalization-group based approach along 



the lines of Ref. 57 can shed new light on this issue. 



Three-nucleon forces (3NFs) are an old but still relevant topic in nuclear 
physics. In spite of many decades of effort, the detailed structure of the 3NF 
is not captured by modern phenomenological 3NF models. Indeed, the global 
analysis presented in Ref. |58| demonstrates that the available models do not al- 
low to significantly reduce the observed discrepancies between the experimental 
data and calculations based on the high-precision NN potentials for breakup and 
polarization observables in elastic nucleon-deuteron scattering. Given the very 
rich spin- momentum structure of the 3NF as compared to the NN force, scarcer 
database and relatively high computational cost, further progress in this fields 
clearly requires input from theory. This provides a strong motivation to study 
the structure of the 3NF within chiral EFT. 

The general structure of the 3NF up to order Q 4 , which also holds at order 
Q 5 , is represented by six topologies shown in Fig. [5] The first non- vanishing 



contributions emerge at N 2 LO (Q 3 ) [59, 60J from the 2vr (a), lvr-contact (d) and 
contact (f) diagrams. The corresponding ttN (ttttN, ttNN, NNN) amplitudes 
at this order are simply given by the Aj = (A, = 1) vertices from the effective 
Lagrangian. The l7r-contact and contact diagrams depend on two LECs cd and 
ce, respectively, whose determination requires few-nucleon data. We will discuss 
the applications of the resulting nuclear Hamiltonian to the properties of few- 
nucleon systems in the next section. It should, however, be emphasized that the 
leading 3NF involves a rather restricted set of isospin-spin-momentum structures, 
which are also included in the phenomenological 3NF models. In particular, the 
longest-range, two-pion exchange topology (a) is well established as one of the 
most important phenomenological 3NF mechanisms. The leading chiral 3NF can, 
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therefore, not be expected to shed new light on the persisting deficiencies in the 
theoretical description of e.g. nucleon-deuteron scattering. The corrections to the 
3NF at N 3 LO generated by the leading loop diagrams are worked out in Refs. 



61 



63 . It is remarkable that the N 3 LO terms do not involve any unknown LECs. 
In addition to the static loop contributions, one also has to take into account the 
1/mjv corrections in the topologies (a) and (d) [63], see also an earlier calculation 



in Ref. 64 . The numerical implementation of the novel N 3 LO contributions in 
few-body studies is non-trivial and requires partial wave decomposition. This 
work is presently in progress, see Ref. [65] for a first step in this direction. 

One of the most interesting features of the 3NF at N 3 LO is, clearly, its rather 
rich spin-momentum structure yielding many operators that have never been 
explored in few-body studies and might potentially be capable of resolving the 
observed discrepancies in 3N scattering. This especially applies to the ring topol- 
ogy (c) and, to a lesser extent, also to the 2ir-lir topology (b). However, the 
observed convergence pattern of the chiral expansion of the 27r-exchange 2N po- 
tential with the leading graphs (b)-(e) in Fig. [2] yielding small contributions and 
the major effect emerging from the subleading diagram (f), see the discussion 
in the previous section, puts a question mark on the convergence of the 3NF at 
N 3 LO. Indeed, since the (large) LECs 02.3,4 saturated by the A-isobar do not 
contribute to the ring and 2tt-1tt exchange 3NF topologies at N 3 LO, one may 
expect that the corresponding potentials of Ref. [62] are not yet converged, see 



also the discussion in Ref. 66 . Thus, one may either need to go to at least N LO 



in the A-less theory or to explicitly take into account the contributions of the 



A-isobar upto N 3 LO. This work is in progress, see Refs. 67 for the first step 
along these lines. 

The parameter-free results (being) obtained in chiral EFT for the various com- 
ponents of the 3NF at large distances rely solely on the spontaneously broken chi- 
ral symmetry of QCD. This opens a very interesting possibility for benchmarking 
with future lattice QCD calculations^] see Refs. 68 69 for first attempts along 
these lines. 

Finally, the four-nucleon force (4NF) also receives its first contribution at N 3 LO 

emerging from tree-level diagrams constructed from the lowest-order vertices of 

dimension A, = 0. The parameter-free expressions for the 4NF at N 3 LO can be 

3 Clearly, care is required to deal with non-uniqueness of the nuclear potentials. The long- 
range part of the 3NF at N 3 LO is, however, uniquely determined after fixing the corresponding 
long-range part of the 2NF. 
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found in Ref. [18]. The contribution of the 4NF to the a-particle binding energy 



was estimated in Ref. 70 to be of the order of a few 100 keV. This provides some 
justification for neglecting four- and more-nucleon forces in nuclear structure 
calculations. 



3 Applications to few-nucleon systems 

Having determined most of the parameters in the nuclear Hamiltonian from 
nucleon-nucleon data, it is now interesting to test it in few-nucleon reactions 
where the ^4-nucleon Schrodinger equation ([!]) can be solved numerically exactly. 
For three particles, the Schrodinger equation is conveniently rewritten in terms of 
the Faddeev integral equations which are usually solved in the partial wave basis, 



see Ref. 71 for details. The Faddeev equations can nowadays be routinely solved 
for any given two- and three-nucleon potentials for both bound and scattering 
states. For a review of recent progress towards including the Coulomb interaction 
in 3N scattering see Ref. (72 j. 

As explained in the previous section, the 3NF at N 2 LO depends on two LECs 
en and ce which need to be determined from few-nucleon data. In Ref. [60] , 
cd and ce were tuned to the triton binding energy and the neutron-deuteron 
{nd) doublet scattering length. The resulting parameter-free nuclear Hamilto- 
nian can then be tested in nucleon-deuteron (Nd) scattering. In Fig. [6j a sample 
of results is shown in comparison with the data. The bands emerge from the 
cutoff variation as discussed in the previous section. Notice that the NLO results 
shown by the light-shaded bands are based solely on the two-nucleon force. We 
further emphasize that the calculations shown do not include the Coulomb inter- 
action and thus correspond to neutron-deuteron scattering. We have corrected 
the proton-deuteron data at E?L = 10 MeV by subtracting out the (estimated) 



Coulomb-force contribution, see Ref. 60 for more details. Remarkably, even 



some rather accurate data for double-polarization observables are available at 
low energy. As a representative example, we show in the top right panel of Fig. [6] 
our results for the nucleon-to-nucleon polarization transfer coefficient K% (N) at 
Eu[ h = 22.7 MeV in comparison with the data of Ref. 73 . More results for dif- 



ferent polarization transfer coefficients at this energy, also based on conventional 



nuclear potentials, can be found in Ref. 74 . For most of the elastic observables, 
one obtains an improved description when going from NLO to N 2 LO which is 
consistent with the description of two-nucleon data at these orders. The increas- 
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Figure 6: Left panel: Differential cross section and tensor analyzing powers T20 
and T21 for elastic Nd scattering at E N h = 10 and 65 MeV. Right panel, top: 
The nucleon-to-nucleon polarization transfer coefficient in elastic Nd scattering 

). Right panel, 



73 



22.7 MeV (proton-deuteron (pd) data from Ref. 
bottom: Nd breakup cross section in the space star configuration (upper data nd, 
lower data pd) . Light (blue) and dark (red) shaded bands show the results from 
the chiral EFT at NLO and N 2 LO, in order. The precise kinematical description 
and references to data can be found in Ref. 1711. 



ing theoretical uncertainty, however, limits the applicability of the N 2 LO chiral 
forces to energies below EF, ~ 100 MeV, where the modern phenomenological 
2NFs and 3NFs also provide an accurate description of the data. There is one 
well-known exception from the generally good agreement between the theory and 
the data at very low energies known as the A y -puzzle. The puzzle refers to the 
strong underprediction of the nucleon vector analyzing power observed at ener- 
gies below E N h ~ 30 MeV for all modern two- and three-nucleon potentials, see 



Ref. 71 for more details. It should, however, be emphasized that A y is (i) very 
small at these energies and (ii) is very sensitive to small contributions to the 
nuclear force |7l| , see Ref. |58| for an extensive discussion. It is, therefore, not 
surprising that the solution to the Ay-puzzle in chiral EFT is not achieved at 
N 2 LO and requires the inclusion of higher-order terms in the Hamiltonian. 

The kinematically very rich deuteron breakup reactions provide even more de- 
tailed insights into nuclear dynamics. At low energies, only very few selected 
observables, mainly the cross section, are available in certain regions of the phase 
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space. While a good agreement between the data and the calculations based 
on the conventional potentials and chiral EFT is observed for the final-state- 
interaction and quasi-free-scattering configurations, large discrepancies occur in 
the case of the space-star configuration (the plane in the CMS spanned by the 
outgoing nucleons is perpendicular to the beam axis, and the angle between the 
nucleons are 120°). This is visualized in the right bottom panel of Fig. [6j where 
a comparison should be taken with the upper sets of neutron-deuteron (nd) data. 
It is remarkable that the existing 3NFs have almost no effect in this observable. 



For recent studies of related breakup configurations see Refs. 75,76 . At higher 
energies the situation is similar to the elastic channels with the predictions from 
chiral EFT being, generally, in agreement with the data but showing a rapidly 



increasing theoretical uncertainty. We refer the readers to Refs. 77-79 , where 
the high-precision cross section and analyzing powers measured recently at KVI 
at Eg h = 65 MeV and covering a large part of the available phase space are con- 
fronted with the theoretical calculations. For a detailed review on 3N scattering 
at intermediate energies the reader is referred to [58| . 

The four-nucleon (4N) continuum provides another interesting and, given the 
appearance of low-energy resonance structures, also very sensitive testing ground 
for nuclear dynamics. It also offers the possibility to probe isospin channels 
not accessible in nucleon-deuteron scattering. The solution of the Schrodinger 
equation for 4N scattering states still represents a major challenge so that only a 
restricted set of calculations, typically at low energies, is available. Interestingly, 
the A y puzzle persists in the 4N system where it becomes even more striking due 
to a much larger magnitude of A y . The very recent study by the Pisa group 80 
shows that, differently to the 3N system, the ^-puzzle in the 4N system is 
significantly reduced by the chiral 3NF at N 2 LO with the LECs cjj and ce being 
adjusted to the 3 H and 4 He binding energies. Given the lack of space, we refrain 
from a more detailed discussion of 4N scattering and refer the reader to the review 



article 80 



The nuclear Hamiltonian at N 2 LO was also applied to compute the spectra of 
light nuclei. In particular, one obtains 7.7 . . . 8.5 MeV and 24.4 . . . 28.8 MeV for 
the triton and a-particle binding energies (BE) at NLO. These results agree well 
with the experimental values of 8.482 MeV and 28.30 MeV, respectively. With 
the triton BE being used as input in the determination of C£> e, the a-particle 
BE at N 2 LO, 27.8 . . . 28.6 MeV, is improved as compared to NLO. 

In all applications discussed so far the LECs entering the 3NF were determined 
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from the triton BE and the neutron-deuteron scattering length. Given the strong 
correlation between these observables known as the Philips line and caused by the 
large S-wave scattering lengths in the two-nucleon system, the resulting values for 
cd,e suffer from a sizable uncertainty. Other possibilities to determine these LECs 
considered in the literature include fitting to the triton and a-particle BEs [81] or 
the properties of light nuclei [82]. Recently, Gazit et al. [83] exploited the fact that 
the LEC cd does not only contribute to the 3NF at N 2 LO but also governs the 
strength of the dominant short-range axial vector exchange current and, therefore, 



can be determined from weak processes. Using N 3 LO 2NF of Ref. [43] combined 
with N 2 LO 3NF with cd,e being calibrated to the triton BE and half-life time, 
the authors of Ref. [83] obtain the a-particle BE of 28.50(2) MeV. Also, the 
point-proton radii of 3 H, 3 He and 4 He found to be 1.605(5) fm, 1.786(5) fm 
and 1.461(2) fm are in an excellent agreement with corresponding experimental 
values of 1.60 fm, 1.77 fm and 1.467(13) fm, respectively. These results provide 
an important and highly nontrivial consistency check of the chiral EFT approach 
by bridging the strong and axial few-nucleon processes. 

Last but not least, the chiral 3NF at N 2 LO is also being extensively explored 



in connection with the spectra of light and medium-mass nuclei, see 84 and 
Ref. [85] for a recent review article describing state-of-the-art calculations within 
the No Core Shell Model, the limit of neutron-rich nuclei and, in particular, 
oxygen isotopes |86| , the properties of nuclear matter and constraints on neutron 
star radii 87 and the puzzle of the anomalously long beta decay lifetime of 

14 C 



4 Nuclear lattice simulations 

A novel scheme to tackle the nuclear A-body problem that combines chiral EFT 
for nuclear forces with Monte Carlo methods that are so successfully used in lat- 
tice QCD or other fields of physics are the so-called "nuclear lattice simulations" , 
also named nuclear lattice EFT (NLEFT). In what follows, we will give a brief 
outline of this approach and present some first results (for a review with many 
references to earlier related work, see Ref. [89]). 
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Figure 7: Schematic drawing of the space-time lattice. The minimal length is 
the lattice spacing a. The side length L in any spatial direction is an integer 
multiple of a. Protons (p) and neutrons (re) reside on the lattice sites. 



4.1 Formalism 

In NLEFT, space-time is discretized in Euclidean time on a torus of volume 
L s x L s x L s x Lt, with L s (Lt) the side length in spatial (temporal) direction, 
as depicted in Fig. [7J The minimal distance on the lattice, the so-called lattice 
spacing, is a (at) in space (time). This entails a maximum momentum on the 
lattice, p m ax = 7T / a , which serves as an UV regulator of the theory. In contrast to 
lattice QCD, the continuum limit a — > is not taken, as we are dealing with an 
EFT and do not wish to resolve the structure of the individual nucleons. These 
are treated as point-like particles residing on the lattice sites, whereas the nuclear 
interactions (pion exchanges and contact terms) are represented as insertions on 
the nucleon world lines using standard auxiliary field representations. The nuclear 
forces have an approximate spin-isospin SU(4) symmetry (Wigner symmetry) [90] 
that is of fundamental importance in suppressing the malicious sign oscillations 
that plague any Monte Carlo (MC) simulation of strongly interacting fermion 
systems at finite density (for a modern look at this symmetry, see Ref. |91|). 
The derivation of inequalities for binding energies of light nuclei in the Wigner 
symmetry limit is given in Ref. [92]. Because of this approximate symmetry, 
nuclear lattice simulations provide access to a large part of the phase diagram 
of QCD, see Fig. [8j whereas calculations using lattice QCD are limited to finite 
temperatures and small densities (baryon chemical potential). Here, we will 
concentrate on the calculation of the ground state properties and excited states 
of atomic nuclei with atomic number A < 12. 
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Figure 8: Nuclear phase diagram as accessible by lattice QCD (yellow area) and 
by nuclear lattice EFT (blue-grey area). On the abscissa, the nuclear density p 
(with pn the density of nuclear matter) and on the ordinate the temperature T 
is displayed. Figure courtesy of Dean Lee. 



We simulate the interactions of nucleons using the MC transfer matrix pro- 
jection method [5j. Each nucleon evolves as if a single particle in a fluctuating 
background of pion and auxiliary fields, the latter representing the multi-nucleon 
contact interactions (for detailed definitions at LO in the chiral expansion, see 
Ref. [93]). We also perform Gaussian smearing of the LO contact interactions 
which is required by the too strong binding of four nucleons on one lattice site. 
More precisely, in a LO calculation using the two independent four-nucleon con- 
tact operators without derivatives ~ (A^iV) 2 , one finds that in the 4 He system 
the ground state is severely overbound and consists almost entirely of the quan- 
tum state with all four nucleons occupying the same lattice site. This is in part 
due to a combinatorial enhancement of the contact interactions when more than 
two nucleons occupy the same lattice site. This effect is partly overcome by 
higher order four-nucleon operators, but it is most efficiently dealt with by a 
Gaussian smearing procedure, which turns the point-like vertex into an extended 
structure. A detailed discussion of this issue can be found in Ref. |93|. Re- 
markably, the aforementioned configurations lead to a new interpretation of the 



phenomenon of clustering in nuclei, see Ref. 94 



Let us come back to the simulation method. To leading order, we start with 
a Slater determinant of single-nucleon standing waves in a periodic cube for Z 
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protons and N neutrons (with Z + N = A). We use the SU(4) symmetric 
approximation of the LO interaction as an approximate inexpensive filter for the 
first to time steps - this suppresses dramatically the sign oscillations. Then we 
switch on the full LO interaction and calculate the ground state energy and other 
properties from the correlation function 

Z A (t) = {* A \exp(-tH)\9 A ) , (24) 

letting the Euclidean time t go to infinity. Here ^> A is the Slater-type initial wave 
function and H is the nuclear Hamiltonian, expressed in terms of the lattice 
variables and lattice fields. Higher order contributions, the Coulomb repulsion 
between protons and other isospin-breaking effects (due to the light quark mass 
difference) are computed as perturbative corrections to the LO transfer matrix. 
This is symbolically depicted in Fig. [9j The perturbative treatment of all these 
effects is justified as for typical lattice spacings a ~ 2 ma, the maximal momentum 
is Pmax — 300 MeV. Note, however, that due to the Gaussian smearing of the LO 
contact interactions, part of the higher order corrections are also treated non- 
perturbatively. If one is interested in the expectation value of any operator O, 



Eq. (24) has to be generalized as 

Z% = {V A \ eM-tH/2) O exp(-tff/2) \V A ) , (25) 

and the ground state expectation value is obtained as the Euclidean time tends 
to infinity (see also Fig. [9J. 

Excited states are calculated from a multi-channel projection MC method, 
briefly described in Ref. [95] . As a first step, we use various improvements in 
our LO lattice action. This is necessary to suppress as much as possible the 
artefacts from the lattice, such that rotational symmetry SO(3) is broken to 
cubic symmetry SO(3,Z) and further artefacts due to the finite lattice spacing 
a arise. To minimize the effect of these, one performs 0(a 4 ) improvements for 
the nucleon kinetic energy and the Gaussian smearing factors of the contact 
interactions. Moreover, all lattice operators at 0(Q 3 ) are included, in particular 
also the ones related to the breaking of rotational symmetry. Their strengths 
can be tuned to eliminate unphysical partial wave mixing like e.g. between the 
3 Si — 3 .Di and the 3 L>3 partial waves. One then uses a set of Slater determinants of 
a large number of different single standing nucleon waves. E.g. for the calculation 
of the spectrum of 12 C, 24 initial standing waves were used and from these three 
states with total momentum zero and J 7 = mod 4 and one state with J z = 2 
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Figure 9: Schematic diagram for the transfer matrix calculation. For the first 
to time steps, the SU(4) symmetric part of the LO action is employed. This 
serves an inexpensive filter to suppress the sign oscillations. Only then the full 
LO action is used. All higher order corrections are included perturbatively as 
shown. The initial and final wave function is a Slater determinant of Z protons 
and N neutrons. 



mod 4 were constructed. Therefore, the correlation function turns into a matrix, 

Z i i(t) = (¥ A \eM-tH)\K) • (26) 

Diagonalization of this matrix of a given ensemble of states with the required 
quantum numbers leads to a tower of states. In such a way, one is able to 
reconstruct the excitation spectrum of any given nucleus. However, due to the 
required computing resources, the ground state and only a few excited states for 
some nuclei have been computed so far. 

The very low memory and trivially parallel structure of the lattice MC codes 
allows to perform simulations, that scale ideally with several thousands of pro- 
cessors. The computational time scales with the number of nucleons A as A 1 ' 7 at 
fixed volume V and with V 1 ' 5 for fixed A. The average sign - which is a measure 
of the severity of the sign oscillations - scales approximately as exp(— 0.1A). Tak- 
ing the calculation of 12 C as a benchmark, the required CPU time for a nucleus 
with spin S and isospin 7 can be estimated as 

, ( - i 41 exp[0.1(>l - 12) + 3(5 mod 2) + 47] , (27) 



X 



CPU 



X CPV x 



and the memory requirements to store the generated configurations are 

2 

orage j ^ } 
-12(7 



X*' ~ X*':?*- { ^ ) cxp[0.1(,l - 1.2) + 3(,S' mod 2) + 17] . (28) 



Therefore, combining high performance computing with the forces derived from 
chiral effective field theory and fixing the parameters in few-nucleon systems 



Chiral dynamics of few- and many-nucleon systems 27 

allows for true ab initio calculations of atomic nuclei and their structure, with a 
quantifiable uncertainty of any observable under investigation. Before presenting 
first result based on nuclear lattice EFT, it is important to stress the differences 
to other ab initio methods. One distinction is that in NLEFT all systematic 
errors are introduced up front when defining the low-energy EFT. This eliminates 
unknown approximation errors related to specific calculational tools, physical 
systems or observables. By including higher-order interactions one can expect 
systematic improvement for all low-energy observables. Another difference is that 
many different phenomena can be studied using the same lattice action. Once 
the action is determined, it can be used to calculate bound nuclei, the ground 
state of neutron matter or thermodynamic properties at non-zero temperature. 
In addition, NLEFT is utilizing several efficient lattice methods developed for 
lattice QCD and condensed matter simulations, including Markov Chain MC 



techniques, auxiliary fields 96,97, pseudofermion methods |98| and non-local 
updating schemes such as hybrid MC [99 - 101| . 



4.2 Results 

So far, calculations in NLEFT have been performed up to next-to-next-to-leading 
order in the chiral expansion of the nuclear potential. At this order, the two- as 
well as the leading three-body forces are present. Consider first the two-nucleon 
system. At NLO, we have nine parameters that are determined from a fit to the 
S- and P-waves in np scattering. Two further isospin-breaking parameters are 



determined from the pp and nn scattering lengths 1 102 . As shown in Fig. 10 



up to center-of-mass momenta of the order of the pion mass, the empirical phase 
shifts in the neutron-proton system are well described and, furthermore, the 
NLO corrections are small. To arrive at these results, a novel method to extract 
phase shifts from finite volume simulations had to be developed (as the standard 



Luscher scheme is not well suited in case of strong partial wave mixing) [104 



Note that due to the perturbative treatment of the higher-order effects, there are 
no contributions to the 2NF at N 2 LO. The three- nucleon force features only two 



LECs, see the discussion in section 2.3, that are determined in our simulations 



from the triton binding energy in combination either with low-energy neutron- 
deuteron scattering in the doublet channel [103 or the a-particle BE [95] . 



The first non-trivial predictions are i) the energy dependence of the pp 1 Sq par- 
tial wave, which agrees with the Nijmegen partial wave analysis up to momenta 
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Figure 10: S-wave phase shifts in the two-nucleon system. From top to bottom: 
Si (np), 1 5o(np), 1 5o(pp). The open/filled symbols represent the results at 
LO/NLO (including all isospin- breaking effects), respectively. The curves are the 
result of the Nijmegen PWA. While the np phases are fitted, for the pp phase 
only the proton-proton scattering length is input. 



of about the pion mass, see the lowest curve and triangle symbols in Fig. 10, and 
ii) the binding energy difference between the triton ( 3 H) and 3 He, 

£( 3 He) - £( 3 H) = 0.78(5) MeV , (29) 

in good agreement with the experimental value of 0.76 MeV [105 106 . The 



theoretical uncertainty mostly arises from the infinite-volume extrapolation. To 
arrive at the result for the tri-nucleon ground state energies, the finite volume 
expression 



E(L) = —BE — — exp(-bL) 

1j 



(30) 



where BE denotes the (positive) binding energy, E{L) the measured energy in the 



finite volume L 3 and a, b are fit parameters, has been utilized (see also Ref. 107 ). 
The ground state energies of nuclei with A = 4,6, 12 where calculated in Refs. 

showing that at N 2 LO one can achieve a precision of a few percent. 



105 



106 



Refining the underlying action as described above and utilizing the multi-channel 
projection MC method, the spectrum of 12 C was worked out in 95 . In Fig. 



11 



we 



show the clean signals of the first few excited states on top of the + ground state. 
Note that due to lattice artefacts the first excited 2 + state is indeed split into two 
states. This will eventually be overcome by choosing a larger basis of initial states. 
In addition to the ground state and the excited spin-2 state, the calculation gives a 
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Figure 11: Extraction of the excited states of 12 C from the Euclidean time 
dependence of the projection amplitude at LO. The slope of the logarithm of 
Z(t)/Z Q +(t) at large t determines the energy relative to the ground state. 

resonance with angular momentum zero and positive parity at —85(3) MeV, very 
close to the 4 He+ 8 Be threshold at —86(2) MeV. Experimentally this threshold 
is located at —84.80 MeV. This first + excitation is the so-called Hoyle state. 
It plays a crucial role in the helium burning of stars heavier than our sun and 
in the production of carbon and other elements necessary for life. This excited 



state of i2 C was postulated by Hoyle 1 1081 as a necessary ingredient for the fusion 



of three a-particles to produce a sufficient amount of carbon and other elements 
needed for life at stellar temperatures. For this reason, the Hoyle state plays also 
a very important role in the context of the anthropic principle, although such 
considerations did not play any role when this state was predicted |109| . The 
Hoyle state has been an enigma for nuclear structure theory since decades, even 
the most successful Greens function MC methods based on realistic two- and 



three- nucleon forces 110 or the no-core-shell- model employing modern (chiral or 
Vi ow k) interactions [82j[84| have not been able to describe this state. In Table [T] 
we show results for the ground state and the low-lying excited states of 12 C at LO, 
NLO with isospin-breaking and electromagnetic corrections included, and N 2 LO. 
For comparison we list the experimentally observed energies. The error bars in 
Table Q] are one standard deviation estimates which include both Monte Carlo 
statistical errors and uncertainties due to extrapolation at large Euclidean time. 
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LO [0(Q°)] 


-110(2) 


-94(2) 


-92(2) 


-89(2) 


NLO [0(Q 2 )} 


-85(3) 


-74(3) 


-80(3) 


-78(3) 


N 2 LO [0(Q 3 )] 


-91(3) 


-85(3) 


-88(3) 


-90(4) 


Experiment 


-92.16 


-84.51 


-87.72 



Table 1: Lattice results for the ground state Of and the low- lying excited states of 
12 C. For comparison the experimentally observed energies are shown. All energies 
are in units of MeV. 
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Figure 12: Results for the 12 C spectrum (Th) and comparison with experiment 
(Exp). The results of the NLEFT at N 2 LO are shown for the ground state, the 
Hoyle state and the lowest-lying spin-2 state. In the simulations, the 2 + state 
is split into the J z = and J z = 2 projections. Only the J z = component is 
shown here. 



Systematic errors due to omitted higher-order interactions can be estimated from 
the size of corrections from O(Q ) to 0(Q 2 ) and from 0(Q 2 ) to 0(Q 3 ). As seen 
in Table [l] and summarized in Fig. |l~2~l the N 2 LO results for the Hoyle state and 
spin-2 state are in agreement with the experimental values. While the ground 
state and spin-2 state have been calculated in other studies, these results are 
the first ab initio calculations of the Hoyle state with an energy close to the 
phenomenologically important 8 Be-alpha threshold. It is important to note the 
energy level crossing involving the Hoyle state and the spin-2 state. The Hoyle 
state is lower in energy at LO but higher at NLO. One of the main characteristics 
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of the NLO interactions is to increase the repulsion between nucleons at short 
distances. This has the effect of decreasing the binding strength of the spinless 
states relative to higher-spin states. We note the 25 MeV reduction in the ground 
state binding energy and 20 MeV reduction for the Hoyle state while less than 
half as much binding correction for the spin-2 state. This degree of freedom in the 
energy spectrum suggests that at least some fine-tuning of parameters is needed 
to set the Hoyle state energy near the 8 Be+ 4 He threshold. It would be very 
interesting to understand which fundamental parameters in nature control this 
fine-tuning. At the most fundamental level there are only a few such parameters, 
one of the most interesting being the masses of the up and down quarks. Such 
investigations have already been performed to unravel the quark mass dependence 
of the deuteron binding energy and of the S-wave nucleon-nucleon scattering 



lengths 111 112 . The impact on the primordial abundances of light elements 
created by a variation of the quark masses at the time of Big Bang nucleosynthesis 



was also studied in Ref. 113 



4.3 Neutron matter 

Matter purely made of neutrons is not only interesting by itself but also of astro- 
physical relevance as various forms of this strongly interacting quantum many- 
body state are realized in the different layers of neutron stars. Due to the Pauli- 
principle, three-body forces are suppressed in neutron matter and, of course, 
there is also no Coulomb repulsion. As for nuclei, one can perform simulations 
for a fixed number of neutrons in a given volume, thus varying the Fermi mo- 
mentum kp (density p) of the neutron matter. For ./V (spin-saturated) neutrons 
in a periodic cube of side-length L, the Fermi momentum (density) is 

(3vr 2 iV)V3 k 3 
k F = , p=— 2 . (31) 

Varying L from 4 to 7 and N = 8, 12, 16 corresponds to a Fermi momentum 
between 88 and 155 MeV, i.e. densities between 2 and 10% of normal nuclear 
matter density, pjy = 0.17 fm 3 . Interestingly, neutron matter at k F ~ 80 MeV is 
close to the so-called unitary limit, where the S- wave scattering length is infinite 
and the range of the interaction is negligible. At lower densities corrections due 
to the (finite) scattering length become more important, and at higher densities 
corrections due to the effective range and other effects become important. In 
the unitary limit the ground state has no dimensionful parameters other than 
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Figure 13: Results for ground state of strongly interacting neutron matter 
Eq/Eq™ 3 versus the Fermi momentum kp using NLEFT at LO (open trian- 
gles) and NLO (filled triangles). For co mpar ison, the r esult s of FP 1981 11201 , 
APR 1998 Fm], CMPR v6 and u8' 2003 [l22], SP 2005 [l23] and GC 2007 ^124 
based on different many-body techniques are also displayed. 



the particle density. Thus, the ground state energy of the system should obey 
the simple relation Eq = $,Eq b for some dimensionless constant £, with E^ 
the energy of non-interacting particles. The universal nature of the unitary limit 
endows it relevance to several areas of physics, and in atomic physics the unitarity 
limit has been studied extensively with ultracold 6 Li and 40 K atoms using a 
magnetic-field Feshbach resonance (a summary of recent determinations of £ can 
be found in Ref. |H4|). There have been numerous analytic calculations of £, 



employing the full arsenal of available many-body techniques, see Ref. 115 for a 



recent review. A benchmark calculation for the four-particle system was reported 
. Around the unitary limit, the ratio Eq/Eq cc can be parameterized 



in Ref. 

as 



114 



Eq 

rpfree 
^0 



s c l 



+ 6 k F r nn + . . . 



(32) 



in terms of the neutron- neutron scattering length a nn and effective range r nn . In 
what follows, we use £ = 0.31(1) and £i = 0.81(1) as determined from lattice 



simulations for two-component fermionic systems in Ref. 116, 117 



In the framework of NLEFT, neutron matter was studied in Refs. |TT8l[Tlg 



The resulting energy of the neutron quantum state compared to a free en- 
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semble is shown in Fig. 13, in comparison to earlier calculations using different 
many-body techniques. For most cases in the range of densities considered, the 



agreement is good. The universal parameter £2 from Eq. (32) is determined to 
be in the range from 0.14 to 0.27 119| . In principle, it can be measured in any 
two-component fermionic system. 

In the future, nuclear lattice simulations of neutron matter can be used to 
investigate the interesting problem of a possible P-wave pairing, as we can dial 
the strength of the neutron-neutron interactions in the various partial waves by 
varying the strength of the corresponding LECs. Obviously, this is not possible 
in nature. Also, simulations for more neutrons and different lattice spacings 
are necessary to get a better grip on the neutron equations of state (EoS) at 
higher densities and also get a better control of the lattice errors. With the 
observation of a neutron star of two solar masses |125| , there exist now much 
more stringent constraints on the neutron EoS (see e.g. Refs. 87, 126) ) and 
therefore such investigations based on NLEFT will be an important tool. 
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